Breakdown of the interlayer coherence in twisted bilayer graphene 
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Coherent motion of the electrons in the Bloch states is one of the fundamental concepts of the 
charge conduction in solid state physics. In layered materials, however, such a condition often breaks 
down for the interlayer conduction, when the interlayer coupling is significantly reduced by e.g. large 
interlayer separation. We report that complete suppression of coherent conduction is realized even 
in an atomic length scale of layer separation in twisted bilayer graphene. The interlayer resistivity 
of twisted bilayer graphene is much higher than the c-axis resistivity of Bernal-stacked graphite, 
and exhibits strong dependence on temperature as well as on external electric fields. These results 
suggest that the graphene layers are significantly decoupled by rotation and incoherent conduction 
is a main transport channel between the layers of twisted bilayer graphene. 

PACS numbers: 71.20.-b, 71.20.Ps, 71.18,+y, 72.20.My 



In many layered systems, the interlayer coupling is 
one of the key parameters for altering their electronic 
properties [l|-l6|- When a thick insulating block is in- 
serted between the metallic layers, the interlayer cou- 
pling can be significantly reduced, leading to breakdown 
of the interlayer coherence, as nicely demonstrated in 
two dimensional electron gas (2DEG) in semiconduc- 
tor superlatticeQ. Such an "interlayer version" of the 
Mott-Ioffe-Regel limit is realized when the layer sep- 
aration exceeds the mean free path across the layers, 
which is evidenced by qualitatively different temperature 
dependences of the intralayer (metallic) and the inter- 
layer (semiconducting) resistivities. The intriguing in- 
terlayer conduction has been observed in various sys- 
tems including high-T c cuprates [l[, organic crystals 
dichalcogenides Q , graphite 0, [B[ and semiconductor su- 
perlattices @, but its underlying mechanism related to 
the interlayer incoherence has been under debate in last 
decades Q. 

Graphene, an ideal 2DEG system, also exhibits rich 
electronic properties depending on how it is stacked on 
top of another graphene layer While bilayer 

graphene in Bernal stacking has massive charge carriers 
with a zero band gap, twisted bilayer graphene with a 
random orientation of the layers has a massless electronic 
dispersion similar to that of monolayer graphene fioj ]. 
Twisted bilayer graphene is of particular interest be- 
cause several intriguing prop erties such as renormaliza- 
tion of Fermi velocity [lO, 11 L van Hove singularities [13], 
and electronic localization [ill. fl3j j were recently discov- 
ered. Experimental studies including angle-resolved pho- 
toemission spectroscopy [l4j . scanning tunneling spec- 



troscopy [l2l . Raman spectroscopy [TBI, [l6| , and in-plane 
transport |l7l. [Hj , suggest that the layers are decoupled 
in twisted bilayer graphene, and the misoriented layers 
are often considered as being electrically isolated. How- 
ever, it is still not clear what sense the layers are decou- 
pled on an atomic length scale of the layer separation, 
and how strongthe interlayer coupling is in twisted bi- 
layer graphene [ll, . 

In this Letter, we present experimental evidence for 
complete suppression of the interlayer coherence between 
in twisted bilayer graphene. We found that the interlayer 
resistivity of twisted bilayer graphene is much higher 
than that of Bernal-stacked graphite by at least 4 orders 
of magnitude. In contrast to the almost temperature- 
independent intralayer resistivity, the interlayer resis- 
tivity exhibits strong negative temperature dependence. 
These results demonstrate that incoherent tunneling is 
a main conduction channel between the layers, due to 
the momentum-mismatch of the 2 dimensional Fermi sur- 
faces from each layer. We thus suggest that the twisted 
bilayer graphene provides a rare example of the layered 
systems showing the interlayer incoherence even with a 
small layer separation of ~ 4 A. 

In order to fabricate the graphene cross junction we 
employed a mechanical transfer process ■ Twisted bi- 
layer graphene is formed in the overlapped region of two 
monolayer graphene layers, as shown in Fig. 1(a). Dur- 
ing the process, the surfaces of both the lower and the 
upper graphene layers are kept free of contact to chemical 
solvents until they overlap with each other. This is cru- 
cial for ensuring tight contact between the layers without 
any chemical adsorbates in-between. 
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the overlapped region (C) with those of the bottom (A) 
and the top (B) layers. The Raman 2D peak at C is blue- 
shifted (2682 cm -1 ) as compared to the peaks at A and B 
(2672 cm -1 ). By contrast, the Raman G peaks appear at 
the same position of ~ 1582 cm -1 for the twisted bilayer 
and the monolayer regions [Fig. 1(e)]. These observations 
confirm that the graphene layers are indeed coupled at 
high energy levels in the overlapped region. According 
to Ref. 0, a blue-shift of the 2D peak by ~ 10 cm -1 
with an enhancement of its intensity corresponds to the 
twist angle of 15°-30°. 

The interlayer transport measurements in the twisted 
layers were done using four branches of monolayers as 
electrical leads. For both bottom and top monolayers, 
the intralayer resistivity as a function of back-gate volt- 
age (V g ) shows a narrow peak at Vg_ = 2 V without sig- 
nificant temperature dependence [21|. For the interlayer 
transport, we apply the current between each branch of 
the bottom and the top layers and measure the interlayer 



FIG. 1: (color online) (a) Optical image of a cross junction 
of two monolayer graphene sheets (A and B) forming twisted 
bilayer graphene in the overlapped region (C). (b) AFM image 
with a scale bar of 500 nm. (c) The height profile along the 
line indicated in (b). (d) The height histogram of the boxed 
region in (b). The data is fitted by the sum (black) of two 
Gaussian functions (red), yielding d ~ 0.4 nm as indicated by 
the arrow, (h) Raman spectra taken from the monolayer (A 
and B) and the twisted bilayer (C) regions. The inset shows 
the Raman map of the relative frequency shift of the G band 
and the 2D band in the twisted bilayer region with respect to 
those of the monolayer region. 



Before discussing the interlayer transport properties of 
twisted bilayer graphene, we confirm that in the over- 
lapped region, two graphene layers are tightly stacked. 
Figure 1 shows the optical and the atomic force mi- 
croscopy (AFM) images of the twisted bilayer graphene 
sample that is mainly discussed in this manuscript. The 
separation between the layers is d ~ 0.4 nm, determined 
by the line-cut height profile [Fig. 1(c)] and the height 
histogram over the step from the bottom to the top lay- 
ers [Fig. 1(d)]. This is very close to the layer spacing of 
d ~ 0.34 nm in Bernal-stacked graphite. 

Twisted stacking of the graphene layers in the over- 
lapped region is also demonstrated by Raman spec- 
troscopy. According to recent studies 0, EH , the double 
resonance peak (2D) of twisted bilayer graphene resem- 
bles a single Lorentzian peak rather than multiple peaks 
as found in Bernal-stacked bilayer graphene. Further- 
more, the Fermi velocity of twisted bilayer graphene is 
known to be reduced as compared to the monolayer in- 
ducing a blue-shift of the 2D peak. Both features are 
found in our sample when we compare the spectra of 
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FIG. 2: (color online) (a) Schematic diagram of a twisted 
bilayer graphene device, (b) Optical image of a twisted bi- 
layer graphene device (SI) with a scale bar of 10 /im. (c) 
Current-voltage characteristics of the interlayer conduction 
at various temperatures from 2 K to 280 K (using the same 
color code in (d)) (d) The interlayer resistivity of twisted bi- 
layer graphene (SI) with variation of the gate voltages (V^) 
at different temperatures, (e) The temperature dependence 
of the normalized interlayer and the intralayer resistivity at 
V g — Vn- (f) Temperature dependence of the pinter(T) for 
several twisted bilayer graphene devices (S1-S4) at V g = Vn- 
The crossover temperature T* from the T-independent to the 
strong T-dependent regimes is indicated by the arrows. 
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voltage drop using the other two branches as illustrated 
in Fig. 2(a). In total, four devices of the twisted bilayer 
graphene were studied (Table 1), and all devices show 
qualitatively similar behaviors of the interlayer conduc- 
tion. We present the results obtained from SI [Fig. 2(b)] 
unless otherwise noted. Figure 2(c) displays the current- 
voltage (I-V) characteristics through the layers at the 
charge neutral point with a back-gate voltage of V 9 =Vn 
in the temperature range from 2 K to 280 K. The I-V 
characteristics exhibit the ohmic behaviors, and it is sus- 
tained down to 2 K. 

Despite the linear I-V characteristics the interlayer re- 
sistivity (pinter) is very large. Figure 2(d) shows the in- 
terlayer resistivity Pinter as a function of V g at various 
temperatures. Here we define pinter = Rinter'A/d, where 
A is the area of the overlapped region, d is the interlayer 
distance, and i?i n tor is the interlayer resistance [13] • At 
280 K, /9inter~ 2000 fkm, which is at least four orders 
of magnitude higher than the reported values of the c- 
axis resistivity p c for graphite. Typically p c ranges from 
~ 10~ 4 57cm for single crystalline graphite [3] to ~10 -1 
ficm for highly oriented pyrolytic graphite (HOPG) 0,0]. 

The interlayer resistivity shows a strong temperature 
dependence as shown in Fig. 2(e). Pi n ter{T) at charge 
neutral point (V g = Vn) increases with lowering tempera- 
ture, and it saturates at ~ 7800 ficm below the crossover 
temperature T* ~ 90 K. The similar temperature de- 
pendence was observed in other samples [Fig. 2(f)]. 
The negative temperature dependence (dpi nter /dT < 
0) is distinct from the intralayer resistance measured 
for the bottom and the top monolayers, which is al- 
most temperature- independent [Fig. 2(e)]. More impor- 
tantly, this behavior contrasts to the metallic behavior, 
dp c /dT > found in single crystalline Bernal-stacked 
graphites Q. These results provide experimental evi- 
dence for the incoherent interlayer conduction between 
the layers 22 1. 

Recent theoretical studies on twisted bilayer graphene 
suggested that Dirac cones from the two monolayers are 
misoriented as shown in Fig. 3(a) [l^]. When the twist 
angle is away from the commensurate angle, the overlap 
of the FS's from each layer does not occur, and the energy 
difference between states with the same extended mo- 



TABLE I: The characteristics of twisted bilayer graphene 
samples. The interlayer resistivity (pinter), the interlayer 
equilibration rate (t^ c ) at 2 K are listed together with the 
crossover temperature (T*) for V g = Vn- The relative shift 
of 2D Raman band (A213) of the twist bilayer graphene with 
respect to that of monolayer graphene is also listed. 
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FIG. 3: (color online) (a) The momentum mismatch of two 
sets of Dirac cones (red and blue) from different layers in 
twisted bilayer graphene. (b) The neighboring Dirac cones 
overlaps at high energies. Thermally-excited tunneling be- 
tween two FS's of each layer is indicated by the arrow, (c) 
Phonon-assisted tunneling with a momentum exchange of q 
enhances the overlap between the FS's. (d) The interlayer 
resistivity (pinter) of twisted bilayer graphene (SI) at dif- 
ferent Vg's as a function of temperature. The solid lines 
are the fit of the metal-insulator-metal junction model whose 
schematic illustrations are shown in the inset (see the text), 
(e) The V g dependence of the conductivity a(0) at zero tem- 
perature limit. The open (solid) symbols represent the case 
of hole (electron) carriers. The height U (f) and the width w 
(g) of an effective barrier estimated from the fitting are shown 
as a function of V g . The solid line is the guide-to-eyes. The 
crossover temperature T* (f) and the resistivity ratio pinter(2 
K)/ Pinter (280 K) (g) are also plotted. The interlayer separa- 
tion d = 0.4 nm is indicated by the dashed line in (g). 



mentum is typically much larger than the Fermi energy. 
In this case, the interlayer tunneling of charge carriers 
with conservation of their in-plane momentum is signifi- 
cantly suppressed. The interlayer equilibriation time e.g. 
the RC time trc for the interlayer conduction is expected 
to be 10~ 6 - 10~ 12 sec for u = 5x 10 12 cm" 2 at zero tem- 
perature limit. Using the measured <r(2 K) at V g = 30 
V and the interlayer dielectric constant eqq ~ 2.6eo [IH, 
we estimate trc ~2x 10 -10 sec, that is well within the 
range of theoretical estimates [19| . 

For quantitative analysis, we employ a model of metal- 
insulator-metal junction [23( 1 . The current through a po- 
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tential barrier is given by 



T(E)D b (E)D t (E + eV)[f b (E) -/*(£ + eV)]dE, 

(1) 

where Db.t(E) and fb,t(E) are the density of states 
(DOS), and the Fermi-Dirac distribution function for the 
bottom and the top monolayers, respectively. In the 
model, the conductivity saturates at a zero temperature 
limit cr(0) - D b (E F )D t (E F )T(0). Assuming no signif- 
icant Independence of T(0), cr(0) ~ \V g — Vm\, which 
roughly agrees with our experiments [Fig. 3 (e)]. The 
finite and relatively large cr(0) near the charge neutral 
point is due to the electron- hole puddles [24}, where the 
charge density inhomogeneity ofn-o^SxlO 11 cm -2 are 
expected [25| , 

We now discuss the temperature dependence of the 
interlayer resistivity, which is determined by tunneling 
probability between the two monolayers. Since the inter- 
layer tunneling occurs when the Fermi circles intersect 
in the extended Brillouin zone [lj|, it is allowed for ex- 
cited carriers at higher energies where the two neighbor- 
ing Dirac cones overlap each other as shown in Fig. 3(b). 
Thus we can consider an effective barrier with its height 
U corresponding to the energies required for interlayer 
tunneling between two misoriented FS's [27j. The tem- 
perature dependence of Pinter(T) is well reproduced by 
the fit using Eq. (1) as shown in Fig. 3(e). Here we as- 
sumed a square barrier with a height of U and a width 
of w [Fig. 3(e)] using the Wentzel-Kramers-Brillouin ap- 
proximation. For the DOS of monolayer graphene, the 
phenomenological DOS is used with a broadening factor 
r = 60 meV gSEilii. 

The resulting barrier height U is reduced at high 
V g — Vn I 's showing electron- hole symmetry, which resem- 
bles the behaviors of T* [Fig. 3(f)]. This is indeed what 
is expected since thermal energy needed for the interlayer 
tunneling decreases when the Fermi circles become larger 
with increase of carrier density. The barrier width w is 
found to be 1.0-1.3 nm with a moderate V g dependence 
similar to that of p in ter{^ K)//9 inter (280 K) [Fig. 3(g)]. 
Note that the estimated w is about 3 times larger than 
the interlayer separation d rs 0.4 nm measured by the 
AFM [Fig. 1]. We tested various cases with differ- 
ent magnitude and shape of U, but w > d 0.4 nm 
is needed in order to explain the observed strong tem- 
perature dependence, i.e. a large value for the ratio 

Pinter 

(2 K) / Pinter (280 K). 

We attribute this discrepancy to the interlayer tunnel- 
ing assisted by phonon at high temperatures [7|. In this 
case, the carriers are scattered by acoustic phonon modes 
within the phonon sphere of diameter q sa ksT/ hv s where 
v s is sound velocity. In particular, in graphene q can be 
much larger than the size of Fermi circle at moderate tem- 
peratures 28(, and thus phonon-assisted tunneling sig- 
nificantly enhances the overlap between the misoriented 



FS's as illustrated in Fig. 3(c). With lowering tempera- 
ture, therefore, not only the thermally-excited tunneling 
but also the phonon-assisted tunneling are suppressed, 
and the effect of the momentum mismatch of the FS's 
become pronounced. This results in the stronger temper- 
ature dependence of pinter {T) , leading to more significant 
suppression of interlayer tunneling than expected in the 
metal-insulator-metal junction model. Further studies 
are highly desirable for understanding the incoherent in- 
terlayer tunneling process assisted by phonon and other 
bosons, as discussed recently @ . Our findings clearly 
manifest that the momentum-mismatch of the FS's due 
to misorientation is essential for breakdown of the inter- 
layer coherence in twisted bilayer graphene. 

In conclusion, using graphene cross junctions we in- 
vestigate the interlayer transport properties of twisted 
bilayer graphene. The large magnitude of the resistiv- 
ity and its strong dependence on temperature and also 
on external electric fields reveal that the main interlayer 
conduction mechanism is incoherent tunneling due to the 
momentum-mismatch of FS's. Thus the twisted bilayer 
graphene is a rare example of the layered systems where 
the interlayer incoherence is induced by misorientation of 
the layers even with an atomic scale of layer separation. 
We envisage that the interlayer conduction can be mod- 
ified by the insertion of other nanosheets using multiple 
stacking [29| or chemical intercalation [3(| ■ Therefore, 
graphene cross junctions with further functionalization 
can provide a new platform incorporating vertical trans- 
port of graphene layers for possible electronic applica- 
tions. 
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